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In this paper it is shown that the Lyman-a transition of a 
single hydrogen-like system driven by a laser exhibits macro- 
scopic dark periods, provided there exists an additional con- 
stant electric field. We describe the photon-counting process 
under the condition that the polarization of the laser coincides 
with the direction of the constant electric field. The theoret- 
ical results are given for the example of ''He'^ . We show that 
the emission behavior depends sensitively on the Lamb shift 
(W. E. Lamb, R. C. Retherford, Phys. Rev. 72, 241 (1947)) 
between the 2si and 2pi energy levels. A possibly realiz- 
able measurement of the mean duration of the dark periods 
should give quantitative information about the above energy 
difference by using the proposed photon-counting process. 

PACS numbers: 42.50.Lc, 42.50.Md, 32.90.-|-a 



I. INTRODUCTION 

For the first time coherence effects in hydrogen-like 
systems were found by the observation of quantum beats 
(Stark beats) in the Lyman-a transition In these 
experiments only the metastable 2s state of atomic hy- 
drogen is initially populated. Switching on a constant 
electric field leads to a build-up of a coherent superpo- 
sition of the upper 2p and 2s levels, and the radiative 
decay shows an interference pattern known as quantum 
beats. 

This uncommon behavior of hydrogen-like systems 
suggests that interesting effects may occur if the Lyman- 
a transition is driven by monochromatic laser light and 
if an additional constant electric field leads to a coher- 
ence between the upper levels 2p and 2s. We show that 
the resulting photon-counting process is similar to the 
one predicted by Dehnielt for a different system with 
two excited states, one rapidly decaying and the other 
metastable, driven by two lasers. Semiclassically one ex- 
pects for the Dehmelt system periods of constant fluores- 
cence intensity due to the strong transition (light period), 
interrupted by periods of zero intensity, while the atomic 
electron is shelved in the metastable state (dark period). 
These photon statistics have been observed experimen- 
tally dl, and the above semiclassical idea has been an- 
alyzed quantum mechanically In an alternative 
experiment Hulet and Wineland proved the existence of 
macroscopic dark periods in the fluorescence intensity of 
a single ion influenced by a magnetic field, when a single 
laser is tuned near one of the principal transition reso- 
nances 0. Hegerfeldt and Plenio Q proposed another 
mechanism for dark periods which is not based on the 



existence of a metastable state. They studied a three- 
level atom with two strong electric dipole transitions to 
one common ground state, driven by a single laser. The 
existence of macroscopic dark periods in the fluorescence 
light is due to a quantum coherence effect. The premise 
is a very small energy separation of the upper levels in 
conjunction with parallel transition dipole moments. Be- 
cause of this an experimental realization of the latter 
physical system seems to be difficult . 

In this paper we discuss the photon-counting process 
of a single hydrogen-like system, driven by a single lin- 
early polarized laser and additionally influenced by a con- 
stant electric field. As for the Dehmelt system, there 
exists a semi classical explanation for the occurrence of 
macroscopic dark periods in the emission process of the 
hydrogen-like system as follows. The strong Lyman-a 
transition {2p Is) is driven by the laser light. Because 
of this we expect a constant fluorescence intensity. The 
constant electric field leads to the possibility that the 
atomic electron makes a transition from the 2p to the 2s 
energy level. In this case we have zero intensity (dark pe- 
riod), because a dipole transition to the Is ground state 
is impossible. But on the other hand, due to the con- 
stant electric field, there exists the possibility that the 
atomic electron gets out of the 2s back into the 2p en- 
ergy level, and the emission process starts again. By the 
quantum mechanical treatment of the problem we show 
that under the assumption that the polarization of the 
laser coincides with the direction of the constant elec- 
tric field the above semiclassical explanation describes 
the photon-emission process qualitatively [ p^ . 

As will be seen later, one can regulate the mean du- 
ration of the dark and light periods almost indepen- 
dently by varying the intensity of the laser beam and the 
strength of the constant electric field. The mean dura- 
tion of the dark periods depends sensitively on the Lamb 
shift |11| between the 2s i and 2pi energy level. A pos- 
sibly realizable measurement of this mean duration, as 
it was done in the experiments |^,|^ in the case of other 
systems, should give quantitative information about the 
above energy difference. 

We show that there exists a correspondence between 
this system and the above mentioned mechanism of 
macroscopic dark periods without a metastable state [B| . 



II. QUANTUM MECHANICAL DESCRIPTION 
OF THE PHOTON-COUNTING PROCESS 

We consider a single hydrogen-like system without hy- 
perfine structure driven by a single linearly polarized 
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laser with electric field F l and additionally influenced by 
a weak constant electric field F. The laser is supposed to 
be tuned near the 2pi \si transition resonance. The 

2 2 

Hamiltonian in dipole form for the atom interacting with 
the quantized radiation field is given by 



kA 



kA " 



+ eD- (Fi(t)+F). 



(1) 



Here Ha is the atomic fine structure Hamiltonian p3[ . 
We assume the Lamb shift to be incorporated in Ha 
10. As in Fig. |i| the relevant atomic eigenstates with 
positive magnetic quantum number are labeled from |1) 
to |5). For every \i) (i = 1,...,5) with positive mag- 
netic quantum number nii, there exists a corresponding 
atomic eigenstate with the same principal quantum num- 
ber, the same total angular momentum quantum number, 
the same parity and the magnetic quantum number —rrii, 
which we denote by | — i). Then the atomic Hamiltonian 
is given by 



Ha 



iii=i 
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(2) 



where coij is the transition frequency between the states 
\i) and Note that \i\ > \ j\ implies ujn > ujji. 

To describe the photon-counting process one needs 
the probability density t^r; [io, i]) for finding ex- 

actly N photons at times ti < ■ ■ ■ < f^r in the interval 
[to,t]. We assume that the initial state of the complete 
system is |i7)(i7| ® p{to), where \ fl) is the vacuum of the 
quantized radiation field, and p{tQ) is the initial atomic 
density operator. The probability density w is of the 
form (l|-|l3l 

w{ti,. ..,tN; [to,t]) = Tr (^S(t,tN) JS{tN,tN-l) X 

X ■■■JS{h,to)p{to)) , (3) 



where S{t, t') and J are atomic superoperators. Using the 
quantum jump approach |^ (which is essentially equiva- 
lent to the Monte Carlo wave function approach and 
to the use of quantum trajectories |Q) Hegerfeldt ||l^] 
has determined the superoperators J and S{t, t') for an 
arbitrary atom. Considering Refs. [p|JT7| one finds that 

\n){(2\ (g) (^Jp(f)/Tr(-)) is the state right after the de- 
tection of a photon, provided we have a measurement by 
absorption, which we assume from now on. There ex- 
ists a nonunitary atomic operator Uc(t,to), the so called 
conditional (reduced) evolution operator, describing the 
time development of an atom under the condition that 
no photon is observed in the time interval [toj^]; a-nd 



S{t,to) is given by ^(t,to)p(io) = Uc{t,to)p{to)Uc{t,to)^ . 
The conditional evolution operator is generated by the 
conditional Hamiltonian He and with Dy = (i|D|j), the 
generalized damping terms 



r. 



ijkl 



and the decay matrix 
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D,j • Dki\uJki\ 



r= ^ijik\i){k\ 



(4) 



(5) 



|i|.|j|.|fc| = l 
UI,I'=I>|3 



one finds following Refs. that He is given by 

He = HA + eD • (Fi (t) + F) - iar. (6) 
For the reset operator J one obtains 

JP^ {^J^kl+rklJ^){^\p\k)\J){l\. (7) 



ijkl 
i|>l3 
fc|>|i 



The above operators He and J are rather complicated 
in the case of our ten level hydrogen-like system. There 
are three main difficulties in calculating the probability 
density Eq. (|^) as follows. The superoperator J gen- 
erally carries pure states into statistical mixtures of the 
two ground states |f), | — 1), with in general non zero 
off-diagonal matrix elements. The state right after the 

detection of a photon |/?)(/?| (g ^ Jp(i)/Tr(-)^ generally 

depends on p{t) The superoperators J and S{t,t') 

act on a f 00 dimensional vector space. In the next section 
we show a way out of these difficulties. 



III. SYMMETRY CONSIDERATIONS 

Here we make the assumption that the laser light 
is linearly polarized with polarization in the same di- 
rection as the constant electric field, F = Fez and 
F^it) = Fl cos{ujLt + Vl)gz- Therefore we have an in- 
variance of the Hamiltonian Eq. (|^) with respect to the 
group 0(2) of those orthogonal transformations which 
leave the z-axis invariant. The results of the following 
symmetry considerations are given at the end of this sec- 
tion. 

By T we denote the standard double-valued represen- 
tation of the group 0{2) on the atomic Hilbert space. 
The invariance of the Hamiltonian Eq. ([^) leads to 



T{g) (jp) T{g)^ ^ J {T{g)png)^) 



(8) 



for every group element g and to the analogous equation 
for S{t,t'). Special elements of 0(2) are given by ^{(p), 
the rotation through the angle if about the z-axis, and by 
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S, the reflection in the y-z plain. One can show that the 
projector P'^"-' onto the subspace of the scalar operators 
with respect to the group 0{2) is given by |^ 

i,j>0 

x(N)01+£.e*|-z)(-j|), (9) 

where p is an arbitrary atomic operator and Si is a phase 
factor with the property T{S)\i) = e^l — i). Using the 
cyclic invariance of the trace and Eq. (||), Eq. we 
obtain 

w[ti, . . . ,tN; [to,*]) = Tr {s[t,tN)JS{tN ,tN-l) 

X ••• J^(ti,io)P^"V(io)) . (10) 

Since the atomic state right after the detection of a pho- 
ton is a statistical mixture of the ground states |1), | — 1) 
one finds 

jp(")p = Tr ( Jp) pr (11) 
with the scalar operator 

P. = ^(|l)(l| + |-l)(-l|)- (12) 

From Eq. (||) we also know that J and S{t,t') leave 
the subspace of the scalar operators invariant, and using 
additionally Eq. (^l]) we come to the conclusion that the 
probability density Eq. (||) factorizes into single-photon 
probabilities 

w{ti, ...,tN; [to, A) = Tr (^S{t,tN)pr^ 

xTr (^JS{tN, tN-l)pr) • • • Tr ( J5(t2, tl)Pr) 

xTY(^JS{h,t„)p(°^p{to)) . (13) 

From this we see that the photon-counting process is 
governed by Po{t), the probability density of counting 
no photons until t starting with the ground state, since 

-Po{t - t') = Tr (^JS{t, t')pr^ 0. The conditional evo- 
lution operator Udt., t') is a scalar operator and from Eq. 
(^ one obtains that Uc{t,t') leaves the atomic subspace 
generated by the states with positive magnetic quantum 
number invariant. Using the definition of pr, the relation 
T{S)\l){l\T{Sy = |-l)(-l|,thefact that Uc{t,t') com- 
mutes with T{S) and the cyclic invariance of the trace 
we obtain 

Po{t) = Tr(^S{t,0)pr) 

= TT{Ua{t,0)prUS,Oy) 

= \mt,omr- (w) 



Therefore we only have to consider the atomic states with 
positive magnetic quantum number. By a similar reason- 
ing it is easy to check that the latter is also true in the 
more general case of the first single photon probability 
Tr(J5(ti,to)PWp(to)) inEq. (|l|). 

Because of our symmetry assumption the original ten- 
level atom behaves like a five-level system with a single 
ground state . In the special case of the spontaneous 
Lyman-a transition influenced by the Stark effect this 
agrees with the previous results in literature (see for ex- 
ample [Q). Starting with the state right after the de- 
tection of a photon even the state |5) drops out, because 
it is coupled neither by the laser nor by the constant 
electric fleld. Going over to an interaction picture the 
explicit time dependence of He in Eq. (^) vanishes. We 
introduce the operator 
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(15) 

in matrix form with respect to the atomic basis |1), |2), 
|3) , |4) , where 7 is the Einstein coefficient of the Lyman-a 
transition, il^ = -1^^(21211) is the real Rabi frequency of 

the laser with respect to the 2ni — > Isi transition, Q = 

22 

■|F(2|z|3) is the analogous real constant of the constant 
electric field and A^ = — loh is the detuning of the 
laser with respect to the state Then we finally have 
Poit) ^ \\e-^'\l)r M- 

IV. THE EMISSION BEHAVIOR 

We assume the laser to be tuned near the 2pi — > Isi 

2 2 

transition resonance, such that IA2I < 7. In this case 
TtjAsI is essentially given by the Lamb shift between the 
2s 1 and 2pi energy level, and IA4I is a fine structure 
frequency, which leads to IA3I <C IA4I. We also as- 
sume the electric fields to be chosen to satisfy the relation 
\Q\ <^ I^IlI < lAaj. We approximate the function Po{t) 
by means of a perturbative approach based on the book 
of Kato 1^^. Since additionally the 2pa level couples 
weakly to the laser, we have two perturbation param- 
eters fl/As and 17l/A4. Putting = in Eq. (|l|) 
we denote the resulting operator by M*-"^ , and we define 
A/(^^ = M — AI^^\ There exists a decomposition of the 
form 

M(°) = Af^P^°^ + ) + Pf ^ , (16) 

where Aj*'"' = — iAs is one of the eigenvalues, another one 
is given by 

Af«|-iA4 (17) 
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in first order in 0.^/1^^ and P'f'\ i = 2, 3 are the respec 



tive eigenprojectors, and we define pf*^ = 1 — P2 ~^ 3 



(0) p(0) 



The operator A^''^ is chosen to commute with pf^\ 



We 



obtain for A^"'' in first order in ^II / A4 the result 



(0) 



?^^(|2)(1| 



>(2|)-i?i( 



iA^ 



|2)(2|. 
(18) 



Note that the right hand side of Eq. (nsl) is the con- 
ditional Hamiltonian of a two- level atoni~p|,p7t . In an 
analogous manner we can decompose the operator M in 
the form 



M = AiPi + A2P2 + A3P3 
with each P,; corresponding to P^'^^ , and we have 



e-^^*|l) =e- 



-e-^^*P2|l) 



c-^^'*P3|l). 



(19) 



(20) 



The main idea of our perturbative approach is to approx- 
imate Ai, A2, A3 and the respective projectors separately 
with the aid of Ref. . First of all we are interested in 
the behavior of Po(t) assuming t 3> 7^^- By using Eq. 
( p^ ) and Eq. (18) one can verify that the first and the 
third term in Eq. ( |20| ) decay exponentially on the time 
scale 7~^ while in first order in fi/Aa the real part of A2 
vanishes. Because of this we only have to approximate 
Re(A2) and |jP2|l)|p. In first order in fi/As we have 



P2II) 



-P 



(0) 



- A 



(0) 



1 - P 



(0) 



(21) 



With the definition of the complex number 

a = 1 - — 
A4 

^^iA2 ^ 

2A2A4 
we obtain 
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2A3 


A4 " 
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2A3A4y 



0202 

-2Rc(A2)t " "l 



l-3s! + 2A. 
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4A5 



4A| 



(23) 



assuming t ^ 7 ^ , up to small relative deviations of the 
order fl/flL. In second order in il/A^ one finds 



Af - 



- A 



(0) 



1 - P 



(0) 



m(i)|3) 
(24) 



and this leads to 



ReAa 
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4A2A3 
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(25) 



It is comparatively easy to describe the behavior of Pq (t) 
on the other time scale t ^ 7~^- In zeroth order in 
ilL/A4, il/As we have 



\1] 



where A^^'' can be approximated by Eq. (Il8| 



(26) 



We introduce a time Pq such that 7 ^ To ^ 
(2ReA2)"^ Then for t < Tq the function Po(t) is 
governed by the behavior of a two-level atom with a 
strong transition, while in a large time interval around 
To it is very small, though not vanishingly small, and 
slowly varying. An interruption of the atomic fluores- 
cence longer than Tq is called a dark period. The above 
results concerning Pq (t) guarantee the occurrence of light 
and dark periods in the resonance fluorescence of the 
atom (see for example [D). Following Refs. |^,^ we can 
calculate the mean durations Pl, To of the light and 
dark periods and the probability p for the occurrence of a 
dark period. One finds p = Po(To) and To = (2ReA2)"\ 
which is given by Eq. (|3|), Eq. ( ^ ) respectively. The 
value of Pl can be obtained from P^ = T^/p, where tl 
is the mean time between two photons in a light period. 
This is intuitively obvious, since p^^ is the mean number 
of photons in a light period. We have 



TL 



4Ai 



1-p 



7 



(27) 



by using Eq. (Eq). Thus one finds 



Td = 



3^ 



4A^ 



372 

4A! 



and 



Tl 



4At|a|'(72 



2nl 



4Ai) 



7mf72 



9 A2 

^A4 



972 

4A^ 



(28) 



(29) 



All the above mean values can be obtained from a single 
trajectory of the photon-counting process. 

As a typical example of the occurrence of macroscopic 
dark periods in the emission process of ''He^ we discuss 
A2 = 0, P = 3.6 X 10^ ^, Fl = 2.9 X 10^^, which means 
n = 0.0257, = 57. One finds 



Td = 1.1 X 10" 



Tl 



4 X lO"'^ s, 



(30) 



and for the mean number of photons in a light period 
we obtain p-i = 2 X 10 . Since the Lyman-a transition 
is remarkably strong with a lifetime of about 0.1 ns, one 
has a high fluorescence intensity in a light period and 
a different time scale in comparison with the Dehmelt 
systems in Ref. H. 
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Under consideration of our premises with respect to the 
electric fields we know that / To is almost independent 
of f2. On the other hand To only depends weakly on fi^, 
which is intuitively obvious. As a conclusion, one can 
regulate the emission process with the aid of the electric 
fields. 



V. DISCUSSION 

From the calculation above we have seen that macro- 
scopic dark periods occur in hydrogen-like systems like 
''He^ provided the external electric fields are suitably 
chosen. One might wonder, however, in which way one 
can reach a dark period. As the quantum mechanical 
calculation shows and the intuitive explanation in the in- 
troduction suggests, in order to reach a macroscopic dark 
period the system must be mostly in the 2s state. 

One might be tempted to argue in a simplified way as 
follows. One could assume that the coherent evolution of 
the atom is started by the absorption of a Is — > 2p pho- 
ton, and terminated by spontaneous emission into this 
channel. In order to evolve to an extended dark period 
spontaneous emission must not occur for many lifetimes. 
If we take the effective Rabi frequency = I/4O7 of 
the Stark field one would estimate that the probability 
to obtain an even mixing of 2p and 2s is smaller than 
exp(— 20), or 2 X 10^^, and at this point the atom is 
not dark at all. From the 0.1 ns lifetime of the 2p state 
of "^He"^ one estimates a lower limit of > 0.1 s for 
the mean time of a light period in the emission process. 
This result much exceeds the previously calculated value 
of = 4 X 10""* s from the quantum mechanical de- 
scription, and macroscopic dark periods should be very 
seldom. 

At this point we have to remember that it is the rela- 
tive weight of the 2s state in the emission-free subensem- 
ble that counts rather than the absolute population. 
There are two mechanisms that make the 2s state become 
rapidly predominant in the emission-free subensemble as 
follows. The relative weight of the 2p and the Is state in 
the emission-free subensemble decreases rapidly on the 
time scale 7~^, because those atoms with a spontaneous 
emission from the strong Lyman-a transition leave the 
emission-free subensemble and do not contribute. On 
the other hand the 2s state is metastable and weakly 
coupled to the 2p state. Therefore if the atom is once in 
the 2s state it stays with a high probability, and it re- 
mains in the emission-free subensemble for a long time. 
As a conclusion we obtain the possibly astonishing re- 
sult that the 2s state becomes predominant fairly quickly 
in the emission-free subensemble although the absolute 
population of this metastable state is very small. An es- 
timation of the population dynamics in the emission-free 
subensemble from a very simple rate equation model is 
given in the appendix. 

For the mechanism of quantum jumps in hydrogen- 



like systems there exists a close relation to the proposal 
of macroscopic dark periods without a metastable state 
by Hegerfeldt and Plenio ||^ as follows. If we neglect 
the existence of the weakly coupled 2p3 level, and if we 
consider the conditional Hamiltonian, which is given in 
matrix form with respect to the atomic orthonormal basis 
|1), ^(|2)-f|3)), ^(|2)-|3))by 



Hr = h 
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2V2 
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A2+A3 
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A3-A2 



2\/2 



-A2 



2 

A2+A3 
2 



(31) 



then Eq. corresponds directly to Eq. (5) of 

H, except that the off-diagonal frequency shift terms 
(A3 — A2) /2 (half of the negative Lamb shift frequency) 
are absent j2^]. As a conclusion, the physical system of 
Hegerfeldt and Plenio behaves like a hydrogen-like sys- 
tem without quantum electrodynamical corrections of the 
atomic spectrum. 

In the case of our realistic hydrogen-like system the 
photon-counting process is governed by the detuning A3, 
where TijAsj is essentially given by the Lamb shift be- 
tween the 2si and the 2pi energy level. We assume a 
possibly realizable measurement of the mean duration of 
the dark periods, as it was done in the experiments 
for Dehmelt systems. In this case we can calculate the 
detuning A3 with the aid of the approximation Eq. ( [2^ ) 
to high accuracy by solving Eq. (|2^) for A3 which leads 
to a polynomial equation of sixth degree. If the other pa- 
rameters in Eq. (|28| ) are known, this provides a detection 
of the Lamb shift by using the proposed photon-counting 
process. 



ACKNOWLEDGMENTS 

I am grateful to G. C. Hegerfeldt, D. G. Sondcrmann, 
and A. W. Vogt for helpful discussions. 



APPENDIX 

In this appendix we show that, in the emission- free 
subensemble, the relative weight of the 2s state becomes 
rapidly predominant on the time scale of the inverse 
Lyman-a Einstein coefficient 7^^. This can be seen 
by a simple rate equation model for the emission-free 
subensemble as follows. For convenience we neglect the 
existence of the weakly coupled 2p3 level. By Pi{t) 
{i = 1, 2, 3) we denote the probability that no photon has 
been detected until t and that the atom is in the state 
|z) at time t. We note that t > implies Yl^i=i ^ii^) ^■ 
By Rg, Rji we denote the transition rates due to stimu- 
lated emission of the blue transition 2p — > Is and the red 
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transition 2s — > 2p respectively. For this subensemble 
one has the rate equations 



Pi = -RbPi+RbP2, 
P2 - RbPi -{i + Rb 
Pz ^ RrP2 - RrPz- 



Rr)P2 + RrPs 



(32) 
(33) 
(34) 



The only difference to the usual rate equations [Q is that 
those atoms with a spontaneous emission from the blue 
transition leave the emission-free subensemble. Therefore 
the term 7P2 in Eq. (|3|) is absent. This leads to a decay 
of Pi{t) in time. Assuming R^ <^ Rb,1 and starting 
with the ground state one easily obtains the perturbative 
expressions 



P2{t) = 



R 



+e 



fj.lt 



(35) 



Ml - M2 

RrRb I e 



Ml - M2 VMi - M3 M2 - 
^13* / RrRb RrRb 



, (M2 - M3)(mi - M2) (mi - M3)(mi - M2) / 
where the real numbers /ii are defined by 



Ml/2 = ^ ( (7 + -Rfl 
M3 = Rr- 



2Rb) ± ^i^ + RRf + AR 



We note that mi > M2 ^ Msi and we see that the pop- 
ulation P3{t) of the 2s state |3) in the emission-free 
subensemble increases rapidly from ^3(0) = on the 
time scale M1/2 ^ ^b^ ^ ^'^'^ then remains on a 
low level for a long time of the order fJ,^^ ^ Pr^- On 
the other hand the population Pi{t), P2{t) of the Is, 
2p state respectively decreases rapidly on the time scale 
Rg^ ~ 7"^ so that P3(t) > P2{t),Pz{t) is reached fairly 
quickly. This behavior can also be seen in Fig. I Be- 
cause of the normalization (X]i^i(^) — weight of the 
emission-free subensemble) the conditional probabilities 
are Pj{t)/ J2i P-ii^) U — li2,3). A similar behavior is 
obtained in the quantum mechanical calculation of the 
paper, but a quantitative agreement is not easily avail- 
able from this simple rate equation model. 
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FIG. 1. Relevant energy levels of "'He"'". The fine structure 
frequency between the 2p3 and the 2pi energy level is given 

by 1.75 X 10^^ Hz, while the Lamb shift frequency between the 
2si and the 2pi energy levels is 1.4 x 10^° Hz. In addition 

7 = lO'*^ s~^ is the Einstein coefficient of the Lyman-a transi- 
tion. Note that the above Lamb shift splitting is appreciably 
larger than 7. 

FIG. 2. Estimation of the expected population dynamics 
by means of the simplified rate equation model in the case of 
the parameters Rb = 57, Rr = O.O57. The dashed line, fat 
solid line, thin solid line indicate the population -Pi(i), P2{t), 
P:i{t) respectively. The time axis is given in natural units of 
the inverse Lyman-a Einstein coefficient. 
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